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Abstract 



We apply the time-dependent local-spin-density approximation as general 
theory to describe ground states and spin-density oscillations in the linear 
response regime of two-dimensional nanostructures of arbitrary shape. For 
this purpose, a frequency analysis of the simulated real-time evolution is per- 
formed. The effect on the response of the recently proposed spin-density waves 
in the ground state of certain parabolic quantum dots is considered. They 
lead to the prediction of a new class of excitations, soft spin-twist modes, with 
energies well below that of the spin dipole oscillation. 
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Recent advances in semiconductor technology nowadays allow the fabrication of nanos- 
tructures with many different shapes. In these systems the electrons, which are laterally 
confined at the semiconductor boundary, form a two-dimensional quantum dot with a shape 
which, to a certain extent, follows that of the nanostructure. This opens up the exciting 
possibility to produce and study an enormous variety of quantum dots, or artificial atoms 
as they are often called. For instance, it has been shown that the electronic structure in 
the small vertical quantum dots of Ref. |IJ is given by the successive filling of shells obeying 
Hund's rules as in atoms. Very relevant information about electronic excitations in quan- 
tum dots is also presently obtained from sophisticated far-infrared absorption g and light 
scattering experiments ||. 

Up to now, the great majority of experimental and theoretical efforts were focussed on 
quantum dots with circular symmetry. Many of the properties of circular dots are well 
reproduced by considering the electrons as confined by a parabolic potential, or by a simple 
jellium disk. To treat the electronic interactions, besides exact diagonalization for very small 
dots [|J, the most succesful approaches have been mean field theories like Hartree-Fock (HF) 
and density functional in the local- spin- density approximation (LSDA) . 

The latter ones have been extended using the random-phase approximation (RPA) to 
analyze collective excitations To our knowledge, all theoretical approaches addressing 
collective excitations in 2d quantum dots are limited from the start by the circular symmetry 
assumption. In this Letter we show how LSDA can describe both ground state and linear 
response of 2d quantum dots of arbitrary shape by using, respectively, energy minimization 
and real time simulation of the spin-density oscillations as basic principles. We will show how 
from the response frequencies in the different channels (density, spin and free responses) it 
is possible to gain information about the system deformation in a quantitative way. Besides, 
we will also analyze the effect on the response of the recently proposed spin-density waves 
in the ground state of particular parabolic quantum dots. Static spin-density waves could 
manifest their existence by means of soft spin-twist modes, at energies well below that of 
dipole spin oscillation. 

Several authors have recently addressed the problem of describing quantum dot ground 
states within LSDA. In particular, in Ref. || the single particle Kohn-Sham equations for 
electrons in a parabolic potential were solved avoiding any symmetry restriction by using a 
plane-wave basis. We use here the same LSDA functional of Ref. ||, i.e., the local functional 
based on the von Barth-Hedin |J interpolation of the Tanatar-Ceperley [|l(| results for the 
non-polarized and fully polarized 2d electron gas. However, we employ a different technique, 
based on the discretization of the xy plane in a grid of uniformly spaced points. For each 
spin (77 =T, j) the Kohn-Sham equations read 



i V 2 + v {conf \v) + v^ H \v) + v[f c \v) 



Viv( r ) (!) 

where t>( con /)(r) and v^ H \r) = f dr'p(r')/|r — r'| are, respectively, the confining and Hartree 
potentials. The exchange-correlation contributions are obtained from the local energy den- 
sity £ xc (p,m) by 

v^(v) = -^£ xc (p,m). (2) 
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We have defined total density and magnetization, in terms of the spin densities p v {t) = 
J2i \(Pir)(r)\ 2 , as p = pi + pi and m = — pj, respectively. 

As a test of the numerical code using the xy grid we have checked that for a circular dot 
(namely, the parabolic one confined by imc^r 2 , with r the radial coordinate, ujq = 0.25 H* 
fTT| and N = 20 electrons) we find the same solution that is obtained by solving only 
the radial equation and imposing e~ lliQ as the angular part of the single particle wave 
functions. Next, we have considered different confining geometries for dots with = 20 
electrons. In particular, we present here results for a deformed parabola \m{uj 2 x x 1 + ujyH 2 ) 
with u y = 0.75 u x = 0.22H*; for a square jellium with r s = 1.51, that corresponds to side 
length L = 11.96a^; and for a rectangular jellium with the same r s and sides L y = 0.75 L x = 
10.37a^. A more systematic investigation of ground states for different sizes and different 
confining geometries is left for future work. Here our aim is mainly to show the feasibility 
of the method and to concentrate on the spin-density oscillations. 

The previously mentioned ground states are shown in Fig. 1. In the deformed parabola, 
the density shows an ellipsoidal shape with an aspect ratio similar to u y /u x . Three rows of 
local maxima can be seen in the inner part of the dot, aligned with the long axis. For the 
square jellium we obtain a rather abrupt electron density, with maxima at the corners and 
four additional inner maxima following the square symmetry. A similar structure is seen 
for the rectangle. Quite interestingly, while for the deformed parabola the magnetization 
vanishes everywhere, for the square and rectangle there is a magnetization wave in the 
ground state. The amplitude of this wave is approximately 15% and 25% of the maximum 
density, respectively. This finding is similar to the spin density waves predicted by Ref. 
in some circular parabolic dots. 

The description of spin-density oscillations in quantum dots has raised great interest, 
mainly due to the manifestation of these modes in far-infrared absorption and in Raman 
scattering experiments We refer here to general spin-density oscillations. When both 

spin components oscillate in phase they produce density modes and when they are out of 
phase, spin modes. In circularly symmetric dots, density modes have been studied using the 
Hartree [12] and Hartree-Fock || methods. More recently the LSD A to density functional 
theory has been used in circular dots to describe density and spin channels ||, taking 
into account the coupling between both. All these methods are based on the perturbative 
treatment of the response by diagonalization of the residual interaction within a space of 
particle-hole excitations. They share as essential ingredient the angular momentum selection 
rules given by the circular symmetry. In fact these methods use the well known RPA, based 
on different ground state theories. 

To describe spin-density oscillations in an arbitrary structure the RPA approach becomes 
practically unfeasible because of the enormous dimension of the matrices. This is due to the 
lack of symmetry, which forces to deal with matrices xi x i Hi x ' ■> y'\ m the the formal RPA 
equation \ — X^VphXi where x^ is the independent particle correlation function and V p h is 
the residual particle-hole interaction. The calculation is also complicated due to the breaking 
of degeneracies with deformation, that greatly increases the number of different particle-hole 
pairs contributing to x ■ An alternative approach that permits to overcome these problems 
is based on real-time methods. These originate in fact from the time-dependent HF theory, 



and have been applied with success to nuclear |13| and to cluster physics [pp . In what 



follows we briefly comment this approach and show how it applyies to 2d nanostructures. 
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In the small amplitude limit, it is well known that both real time (TDHF, TDLSDA) 
and RPA methods based on the corresponding ground states become equivalent. 

We have performed TDLSDA calculations by integrating the time- dependent Kohn-Sham 
equations 

d 

where h v is given by the square bracket in Eq. ([I]). We have integrated Eq. (^j using the 
Crank-Nicholson algorithm (for time step n to n + 1) 

(i + ^r 11 ) vsr 1 ' = (i - ^r>) ^ • « 

This is in an implicit problem for Lp^ n+l \ since h^ +1 ^ depends on the orbitals through the 
density and magnetization. In practice this forces to proceed by iteration: with h v from the 
previous time step, a first guess of the new wave functions is obtained solving ([|). These 
are then used to build a new h v and restart iteration. Using a rather small time step it 
is enough to make a double solution for each time step. If h v were constant in time, the 
algorithm would be exactly unitary. We have found that with small At norm conservation 
is fullnled with excellent accuracy. 

In order to excite the oscillation modes of the system an initial perturbation of the wave 
functions is needed. Physically, this corresponds for instance to the interaction with a short 
laser pulse or with an appropriate projectile. In the calculation, it can be mimicked simply 
by a rigid translation of the wave functions by means of the operator T(a a ) = e~ a<T ' v or by 
an initial impulse with II(q (J ) = e~ 1Cia ' r . When either a CT or q CT are small these perturbations 
induce dipole oscillations predominantly and the system's response is restricted to the linear 
regime. Total density and spin modes are obtained with the rigid translations = = a, 
and a-f = — a^ = a, respectively. With the impulse initial conditions these are: = = q 
(density), q-f = — qj_ = q (spin). After the initial perturbation we keep track of the time 
dependent dipole moments dt- (e • r) t for density modes and (e • rcr z ) t for spin modes. Here 
e corresponds to the direction of the initial perturbation given by a or q. 

A frequency analysis of the dipole signal d t gives the response frequencies of the system. 
Fourier tranform methods can be used for this purpose. However, we have found more 
efficient a method of direct peak fitting to the simulated signal. We perform a least squares 
minimization of x 2 = J2t {dt ~ D(t)) 2 , where D(t) is given by 

N 

D{t) = J2 A n cos(uj n t) + B n sm(u n t) . (5) 

n=l 

In Eq. ([D a fixed number of frequencies is assumed. The minimization yields the amplitudes 
A n , B n and frequencies uj n . Of course, it must be checked that the number iV of frequencies is 
large enough to provide a good reproduction D(t) of the time series d t and that convergence 
with increasing iV has been reached. From the fitted D(t) we obtain D(u) as a discrete set 
of Dirac delta functions. In practice, these are smoothed into Lorentzians and the response 
strength is obtained as S(u) = \D(u)\, or the power spectrum as V{u) = \D(uj)\ 2 . 

We have performed the response calculation in real time for the same dots with 20 
electrons discussed above. An excellent agreement with the RPA calculation was obtained 
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for the circular parabola. Figure 2 shows the three responses for the deformed parabola, 
for which the magnetization m(r) is vanishingly small. The density response has only two 
peaks that coincide with the parabola frequencies uo x = 0.29 and u y = 0.22. This is showing 
that TDLSDA satisfies the generalized Kohn's theorem for a deformed parabola. Since the 
density dipole operator only couples to the center of mass motion, absorption in this channel 
can only take place exactly at the frequencies u x and uj y . The situation is completely different 
for the free and spin responses. They are rather fragmented and lie at lower energies. Peaks 
associated to oscillation along different axis are shown in different line types. The free 
response corresponds to keep the effective confining potential in (|3|), v^ con ^ + v^ H ' + v^ xc \ 
fixed to its static value. It thus models the oscillation of non interacting particles in the 
static mean field. By comparing free, density and spin responses we see the different nature 
of the residual interaction in both channels: weakly attractive in the spin response and 
repulsive in the density one. Figure 2 also shows the simulated time series in each case, as 
well as the fitted signal which on the plot scale superposes to the simulated one. 

Figure 3 shows the corresponding results for the iV = 20 electrons dots in the square 
and rectangular jellium. The density response of the square is still characterized by a very 
dominant peak that, nevertheless, is slightly fragmented. The free and spin responses are 
more fragmented. Quite interestingly, the spin response of the square shows four groups 
of peaks with an almost constant separation of 0.05H*. In the rectangle case, the density 
response is clearly showing that the oscillation frequencies are different in x and y directions. 
The same fact is an additional source of fragmentation for the free and spin channels. 

In circular parabolic dots LSDA predicts spontaneous symmetry breaking ground states, 
of the type of spin density waves, for particular sizes PJI5]. This LSDA spin density wave 
is more pronounced in quasi one-dimensional rings JT^]. In Ref. [ljj similar spin density 
waves, as well as Wigner crystallized ground states, have been very recently reported in 
circular parabolic dots using an unrestricted Hartree-Fock approach. There is, however, a 
continuing discussion about the interpretation and possible relevance of these states P, |T5| , |rT 



We present in the rest of this Letter the TDLSDA result for the linear oscillations of a spin- 
density- wave ground state in a circular dot. 

Finally, we have analyzed the oscillations of a dot with N = 24 electrons in circular 
parabolic confinement (ujq — 0.24 H*), for which there is a static spin-density wave in the 
LSDA ground state [[J. We have found that the dipole oscillations of spin densities are 
quite similar to those obtained from a fully circular model ||. However, the spin density 
wave can sustain a new type of oscillation. It is given by an alternating rotation of both 
spin densities in opposite directions, i.e., a spin twist of the static wave excited with the 
rotation operator 7Z(6 a ), with 6^ = — 0j_ being opposite rotation angles for each spin. The 
frequency of this mode is obtained analyzing the time evolution of the circular currents that 
appear after an initial rotation with TZ(6 a ). Figure 4 shows the strength of the spin-twist 
mode, in comparison with the normal spin dipole mode. Spin-twist modes are very soft, with 
energy well below the spin dipole one. They could signal the existence of static spin density 
waves in circular dots [TS|] . The spin-twist frequency is reflecting the curvature of the energy 
minimum corresponding to the symmetry broken ground state with respect to the circular 
one. We expect that circular systems having strong spin density waves in their ground states 
will also exhibit enhanced spin-twist modes. This happens in circular parabolic dots with 
increasing r s for a fixed iV || (the iV-systematics with fixed r s is less clear) or in quasi one 
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dimensional rings [16 



In conclusion, we have shown that TDLSDA can be applied to obtain the oscillation 
frequencies of nanostructures with arbitrary shape. It leads to the prediction of soft spin- 
twist modes in dots with circular parabolic confining and having static spin-density waves 
in their ground state. 

This work was performed under Grant No. PB95-0492 from CICYT, Spain. 
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FIGURES 



FIG. 1. Two-dimensional surface and contour-line plot of the ground state density (two upper 
rows) and magnetization (two lower rows) of dots with N = 20 electrons in a deformed parabola, a 
square jellium and a rectangular jellium, respectively. The magnetization of the deformed parabola 
vanishes everywhere. See text. 

FIG. 2. Simulated and fitted time series (left), and strength functions in arbitrary units (right) 
for the deformed parabola corresponding to free, density and spin excitations. See text. 

FIG. 3. Strength functions in arbitrary units for a square and rectangular jellium. Free, density 
and spin responses are shown. For the square, only the total spectrum is displayed since x and y 
oscillations are essentially degenerate. 

FIG. 4. First row shows the spin-density wave in a parabolic dot of N = 24 electrons as a 2d 
surface and contour plot. Middle panel shows the spin dipole strength, while lower one shows the 
strength of the spin-twist mode in arbitrary units. 
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